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ABSTRACT 



This is an examination of the open ocean deep water formation process of 
the Greenland Sea and how it is effected by the depth dependent thermal 
expansion coefficient. It is hypothesized that free convection associated with 
parcel instability is possible because of the increase in the thermal expansion 
coefficient with pressure in conjunction with requisite ambient temperature and 
salinity profiles. 

Based on the depth dependence of the thermal expansion coefficient , a 
neutral profile model for parcel instability was formulated and the effects on 
mixed layer dynamics were investigated. In situ profiles for the wintertime 
Greenland Sea were examined for potential parcel instabilities. It was shown 
that only small surface salinity increases associated with freezing could lead to 
deep penetrative convection. Finally, an analysis was performed for regions of 
low stability using climatology from the Levitus data set to determine regions 
most likely for deep convection. 
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I. INTRODUCTION 



There are two main regions of deep water formation, near-continent 
(Antarctic) type and open ocean formation, as discussed in Killworth (1979). 

The Greenland Sea is the site of formation of the largest volume of deep water 
in the North Atlantic (Helland-Hansen and Nansen, 1909). Figures (1) and (2) 
show the Greenland Sea area of interest (Koltermann and Luthje,1989) with the 
significant topographic features and the schematic circulation. 

This paper will examine the open ocean deep water formation process of 
the Greenland Sea. Of particular interest is how open-ocean deep water 
formation is influenced by the depth-dependence of the thermal expansion 
coefficient (qc). It is hypothesized that free convection associated with parcel 
instability is possible because of the increase in the thermal expansion 
coefficient with pressure in conjunction with requisite ambient temperature and 
salinity profiles. The thermal expansion coefficient increases with both pressure 
and temperature. The pressure dependence is more significant for cold polar 
regions than it is for the deep convection of the Mediterranean, as shown by 
Garwood (1991). The dependence of a on pressure, or equivalently depth, is a 
nonlinearity in the equation of state and contributes to a reduction in hydrostatic 
stability (Carmack and Aagaard, 1973; McDougall, 1987). This reduction in 
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Figure 1. Greenland Sea Area Chart Area of interest for conditional 
instabilities. 
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Figure 2. Greenland Sea Topography and Circulation: topographic features 
and schematic of circulation. 
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stability only occurs when the gradient of temperature is negative for z positive 
upwards. 

Parcel instability may occur if a water parcel that is cold and fresh is 
displaced downward into water that is warmer and more saline. This greater 
depth causes increased a for a surface parcel that is displaced downward. For 
a sufficiently large change in a, the parcel density may be greater than the 
density of the surrounding water. This parcel would then be unstable and 
would be accelerated downward by gravity, i.e. a parcel instability would occur. 

The traditional explanation of deep-water formation (Nansen, 1906; 

Mosby, 1959) focused on winter cooling of surface water with the static stability 
reduced by the cyclonic circulation of the Greenland Sea gyre. Killworth (1979) 
hypothesized that with progression of winter cooling the depth of the 
overturning water column increases until top-to-bottom homogeneity is 
achieved. However, there have been no observations to date that clearly 
support vertical homogeneity over the whole water column. 

Carmack and Aagaard (1973) proposed double diffusion as the cause of 
the deep water formation. However, Clarke et al. (1990) showed that the 
double diffusion process can not account for the observed increases in bottom 
water oxygen concentration over expected values for vertical mixing driven by 
double diffusion. 
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Cabbeling (Foster, 1972) and lateral double diffusion have been proposed 
as mechanisms leading to deep water formation, but neither of these processes 
are evident in the chimney data presented by Scott and Killworth (1991) 

(Figure 3). McDougall (1984) showed how diffusion along isopycnal surfaces 
may lead to water mass formation by both cabbeling and thermobaricity. This 
study will not focus on isopycnal diffusion, but on surface-driven buoyancy flux. 

The Greenland Sea stations of interest are located in the center of the 
Greenland Gyre (GG) north of the Jan Mayen Current (JMC). Although 
Killworth (1979) believed that the preconditioning phase does not occur in this 
region, more recently Clarke et al. (1990) pointed out that parts of the Nordic 
Seas are considered more susceptible to deep convection than are other areas. 
The Greenland Sea region circulates cyclonically and is weakly stratified. 

These conditions permit baroclinic instability (Killworth, 1979), and surface 
cooling will induce vertical mixing by creating cooler dense water over less 
dense water. Killworth (1979) hypothesized that the violent mixing initially 
occurs during winter and only in the top layers. When overturning through the 
water column results, the entire water column has characteristics that closely 
resemble those of the deep water. Aagaard and Carmack (1989) found that the 
deep waters in the Greenland Seas are about the freshest of all waters found in 
the arctic and subartic. They proposed that the surface freshwater plays a 
larger role than the intermediate water in the formation process. Gascard 
(1990) stated that deep convection which leads to deep water formation is a 
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Figure 3. Thermistor Chain Data: Thermistor chain data from the Greenland 
Sea representing a chimney event. 
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track distance (3km track marks) 





surface-controlled process involving intermediate and surface waters. The 
short duration of about one week for chimney events suggested by Rudels 
(1989) may explain the paucity of observations of convection that reach the 
bottom in the Greenland Sea. Lazier (1973) showed an isolated example of 
deep convection in the Labrador Sea reaching approximately 1500 m, but it 
disappeared within 2-3 days. 

Including the increase in the thermal expansion coefficient with pressure, 
an entrainment model has been developed by Garwood (1991) that includes 
this nonlinearity in the equation of state. This model predicts greater vertical 
mixing and penetrative convection than do earlier models because of the 
nonlinear buoyancy flux enhancement. In the study here, this process will be 
examined with a conventional one-dimensional model applied to deep oceanic 
mixing. The investigation will show an increase in vertical turbulent kinetic 
energy that may enable deep penetrative convection and the hypothesized 
parcel instabilities leading to formation of deep water. 

Before applying the new mixed layer theory, however, a first step in 
Chapter II was to develop a "neutral parcel" model that predicts a neutrally- 
stable temperature profile assuming a typical two-layer salinity profile. This 
profile is neutral to parcel instabilities, i.e. an isolated parcel of water from the 
mixed layer would have a neutral buoyancy when displaced vertically anywhere 
in the water column. The surface temperature, surface salinity, mixed layer 
depth, and the change in salinity at the bottom of the mixed layer are specified. 
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Then the neutral parcel model diagnoses the change in the temperature at the 
mixed layer bottom and the change of temperature with depth below the mixed 
layer. Thus for a given two-layer salinity profile a temperature profile that is 
neutrally stable for displaced mixed layer parcels is prescribed. 

The next step after development of the neutral parcel model is to study 
the effects of the depth dependent thermal expansion coefficient on mixed layer 
dynamics. A calibrated Kraus and Turner (1967) model with the depth 
dependent thermal expansion coefficient is applied to forced and free 
convection to demonstrate the effect on mixed layer dynamics. 

Chapter III will expand in scope from the mixed layer to the entire water 
column. Parcel instabilities in the water column are examined by moving a 
parcel from the surface or from within the water column over the entire depth of 
the water column. Again, the thermal expansion coefficient is considered to be 
a function of depth and temperatures will be adjusted for heating due to 
compression by the water column. Surface parcels moved in the water column 
will then be tested for sensitivity to changes in initial temperature down to 
freezing and then increases in the initial salinity due to brine injection from ice 
formation. 

Chapter IV explores regional climatology to potential parcel instability, 
determining regions of lower stability. Regional stability is tested by examining 
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the change in buoyancy between a surface parcel and the in situ environment 
at 500 m using climatological temperature and salinity profiles. 



9 



II. NEUTRAL PARCEL STABILITY 



A. DESCRIPTION 



The thermal expansion coefficient, a, is defined as 




( 1 ) 



a *a 0 +ap 



( 2 ) 



where p is the relative density and 0 is the potential temperature. The thermal 
expansion coefficient increases with pressure (p) as a nearly linear function of 
pressure (Garwood, 1991). The accurate calculation of a throughout the 
water column requires knowledge of the ambient salinity and the ambient 
potential temperature versus pressure throughout the water column. The 
surface value of the thermal expansion coefficient is determined by 



The values for the density are found using the International Equation of State 
for seawater (Millero and Poisson, 1981). The neutral-parcel model uses a 
fixed value of a&= 0.01 C. The centered finite differencing scheme is used along 
with double precision due to the small differences in density value. An 
equation similar to (3) is used to calculate alpha at the bottom of the water 




. -p(s 0 / e 0 +Ae,p) +p(5 o ,e o -A0,p) 
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column with an adjustment to the potential temperature for the adiabatic 
temperature increase at the bottom of the water column. 

The adjustment for the adiabatic temperature rise is computed using 
Bryden's (1973) algorithm for adiabatic lapse rate as a function of salinity, 
pressure, and temperature. By prescribing the surface potential temperature, 
the salinity and the pressure at the base of the water column, the increase in 
potential temperature is determined using Newton's convergence technique. 

The temperature increase due to compression is added to the surface potential 
temperature to give the corrected potential temperature at the base of the water 
column. With the calculated temperature and given salinity and pressure at the 
base of the water column, the bottom value of alpha (<x bottofn ) is calculated using 
a form of (3). Equation (3) is modified for <Xt, ottom by replacing the surface 
salinity with the bottom values of salinity, pressure, density and corrected 
potential temperature. Using alpha from the top and the bottom of the water 
column, the value for (a) in (2) is determined using (4). Knowing the constants 

2 _ a(bottom)-a 0 ^ 

of (2), the change in the value of alpha from the top to the bottom of the 
column is specified as a linear function of depth. 

The expansion coefficient for salinity is required for the neutral-parcel 
model to evaluate the change in salinity at the mixed layer interface and in 
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evaluating the temperature values in the modeled "neutral parcel" water column. 
The value of p is determined using (5). The expansion coefficient for salinity, p, 



p.1 6p , 0.799 

p 65 p 



(5) 



is assumed constant throughout the water column. The value of 0.799 psu 1 is 
the average of p values over a 400 bar water column (Gill, 1982). 

With the above information, the neutral-parcel model can now be used to 
prescribe the values for the potential temperature profile as a function of depth 
or pressure. The physical meaning of the neutral-parcel model is that the 
profile is stable everywhere for small vertical displacements, with Brunt-Vaisala 
frequency (N 2 >0). However, a surface parcel that is moved downward will be 
perfectly neutral, remaining wherever it is placed with no restoring force. 

For a neutral buoyancy profile, (6) is used to solve for the potential 
temperature as a function of depth, (7). 



Figure (4) shows profile data for station 82 of the Marginal Ice Zone 
Experiment (MIZEX-87) (Quadfasel et al., 1988). The data from MIZEX-87 was 
collected during the RV "Valdivia” cruise 54 in March-April 1988 in the 



p(2)=p[1 -«W[e-ej-p[5-5j] 



(6) 




(7) 
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On 





Figure 4. Neutral-Parcel Model Run on Station 82: The top two panels are the 
T-S profiles of Station 82. The bottom panels show the neutral-parcel model 
run. 
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Greenland Sea and the Fram strait. One objective of this cruise was to find 
and map an active convective element. Station 82 has representative 
temperature and salinity profiles for the center of the Greenland Gyre. In 
Figure (4), the lower right plot (DIFFERENCE) represents the difference 
between station 82's temperature profile and the neutral parcel temperature 
profile. This "difference" demonstrates how close station 82's temperature 
profile is to a neutral parcel temperature profile. This also demonstrates that 
the center of the Greenland Gyre is close to neutral stability. 

B. EFFECT OF a(p) ON MIXED LAYER DYNAMICS 

To demonstrate, the effect of adding a(p) to a calibrated Kraus and 
Turner (1967) model, three cases are examined: forced convection, free 
convection and a combination of forced and free convection. Using the Kraus 
and Turner model (8) and (9) and the constants listed in Table I determined 




8S_ (F+E-^S 
dt h 



A 

dt 



(8) 



A T'dh 

dT_ Qp dt_ (9) 

dt pCf/l h 
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TABLE 1. CONSTANTS IN THE PARCEL INSTABILITY MODEL 



u* 


1.5 cm/s 


g 


983 cm/s 2 


Ah 0 


1800 cm 


AT ° 


+ 1.52 C 


I “° 


0. 475x10-'* s' 1 


I 3l 


0.265X10' 7 ernes' 1 


Qo 


-0.005 cal/cm 


P_ 


1 g/cm 3 


c p 


1 cal/g/K 
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from the NPS mixed layer model (Garwood, 1977) for arctic conditions, the 
results of forced and free convection for both a constant value of a, and a as 
a function of depth will be solved numerically and compared. In (8), F is the 
rate of water freezing, E is the rate of evaporation and P is the rate of 
precipitation. In (8) and (9), h defines the mixed layer depth, AT is the change 
in temperature at the mixed layer, g is gravity, and 0,,/pCp is the downward 
surface temperature flux. Constants c, and c 2 are dimensionless model 
constants. 

The cooling case (Q 0 ( 0) with continued deepening of the mixed layer 
over time will be considered first. The left hand side of (8) is the rate turbulent 
kinetic energy (TKE) is converted to potential energy (PE) by entrainment (layer 
deepening). The first term on the right hand side of (8) is the generation rate of 
the TKE by wind stirring minus viscous dissipation. The second term is 
damping (Q 0 > 0) of TKE by surface heating or free convection production 
(Q 0 ( 0) of TKE by surface cooling. The left hand side of (9) is the rate of 
change of the temperature at the bottom of mixed layer. The first term on the 
right hand side of (9) is the surface heating and the second term is the 
entrainment heat flux, with a mixed layer temperature decrease(increase) if AT 
is positive(negative). 
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For case 1, forced convection (Q 0 = 0), the wind forcing is the only source 
of turbulent kinetic energy. For this case (8) and (9) simplify to (10) and (11). 



er _ atm 

dt h dt 



dh = erf 

dt aghAT-pgAS 
Note that hAT equals a constant: 



( 11 ) 



AT=T-T[Z=-H) 
8A T = dT dT ( h) 

dt dt dr ) 

f (-«)=0 

dT_ ATdh 
dt" h dt 

h2*I+AT^=0 

dt dt 

4(hA 7)=0 
dr 1 



For a constant over depth and a positive(negative) AT, integrating (11) from h 0 
to h will result in a linearly increasing(decreasing) mixed layer depth as a 
function of time. 
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With the same initial conditions, (11) becomes (12) if a is depth 
dependent. Integrating from h 0 to h results in (13). Figure (5) shows the 
deepening(shailowing) of the mixed layer depth for constant a and depth 
dependent a versions of case 1. The mixed layer deepening(shallowing) is 
slowed over a one day period when a is depth dependent. Furthermore, the 



depth dependent a causes a nonlinear change in depth of the mixed layer as 
noted by the slight curvature of the dotted line for positive AT and the dashed 
line for negative AT in Figure (4). 

Case 2, free convection (u.= 0), the condition of surface cooling (Q 0 ( 0) 
for a constant a and a(z) are compared. Using the Kraus and Turner model on 
the free convection case results in (14). Notice in (14) that a appears in both 



3/?. Ci , 1_ 

3f gh^T^ a 0 -a, 



(12) 



h h ^ f 
gh^T 0 {* 0 +a 



(13) 




(14) 



3f a gh& T 
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FORCED CONVECTION 




uiD ut ipdap 



Figure 5. Case 1: Forced convection mixed layer deepening(shallowing), 
constant a is the dot-dashed(solid) line and a (z) is the dotted(dashed) line. 
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the numerator and the denominator. It is not eliminated because a now has 
different rates at the top and bottom of the mixed layer. The first term on the 
right hand side of (15) is the turbulent heat flux at the surface, TW [. The 
second term on the right hand side of (15) is turbulent heat flux at the bottom of 
the mixed layer, Tw' [. 



When awas assumed constant, the a values in the numerator and 
denominator can be divided out which reduces (14) to 




(15) 




(16) 



And recall from earlier that at z=-h, 



dAT_BT 



dt dt 



Then, 



dAT + ^ jdh _ 
dt dt p C p 




Now integrating from h 0 to h results in 



hA T-hA T 0 =-^—°t 
P c p 



Solving (16) numerically gives solutions in Figure (6). 
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For the depth dependent a case, substitute <x=<x 0 -a^z into (14). The value 
of a in the numerator will be a 0 since it is the surface flux component. The 
denominator a value becomes a 0 + a, h(t) at the bottom of the mixed layer 
(z = -h). Using these substitutions and solving at z = -h yields, 



dt p C p a o +S,/j(0 



) 



Integrating from h 0 to h gives (17). Figure (6) shows that the mixed layer is 



/?A T-h 0 A if. 

pc p J « o +a,/J(0 



(17) 



deepening(shallowing); however, the rate of deepening(shallowing) is much 
slower than in case 1. Note that the deepening(shallowing) is much slower 
than that in case 1. For the a(z) case, the magnitude of reduction in the rate of 
mixed layer deepening(shallowing) is smaller than that in the forced convection 
case. Thus the results show that as in the forced convection case, the slowing 
of the mixed layer deepening(shallowing) is caused by the increase in a with 
time. 
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FREE CONVECTION 




uio in qidap 



Figure 6. Case 2: Free convection mixed layer deepening(shallowing), 
constant a is the dot-dashed(dotted) line and a(z) is the dashed(solid) line. 



22 



time in seconds xlO 4 



Cases 1 and 2 were a simple analysis of how the Kraus an Turner (1967) 
model was effected by simply plugging in a(z). Case 3 will examine how the 



Kraus an Turner (1967) model is effected when a(z) is considered in the 
derivation of (8) and (9). For case 3, the combined effects of both forced and 
free convection (Q 0 ( 0 and u. * 0) are examined. The first term on the right 
hand side of (8), generation rate of wind stirring minus viscous dissipation, can 
be rewritten in terms of turbulent heat flux as shown in (18). The top line of (8) 
can be rewritten as (19) using form shown in (18), the generation rate of TKE 
balance by the vertically integrated buoyant dampening. 



The turbulent heat flux is a function of the depth and is written in the form 
Multiply the turbulent heat flux by ag, substitute in the expressions for the 



surface turbulent heat flux and the turbulent heat flux at the bottom of the 
mixed layer, and integrate over the depth of the mixed layer. These operations 




( 18 ) 




( 19 ) 



rv(z)=rv'(0)(i + ^)-tv / (-^(t) 



(20) 
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result in (21). The left hand side of (21) is the integral average of the buoyancy 
flux over the water column. The first term on the right hand side is the 
buoyancy flux component due to the net downward temperature flux. The 
second term is the buoyancy flux component due to the bottom heat flux. 

*|)(-g)^/.:«9(|)(A rf )* (21 ) 

For constant a, we can reduce (21) by dividing out the gravity and solving 
the integral to give 

(22) 

hJ-h p Cp dt 2 

For the depth dependent a, we will substitute into (21) the form a = a 0 - a,z 
from (2): 



f °agT'w'dz 

J -h 

+ /;s(« 0 -a,^(|)( 

Expanding the integral gives (23). 
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(23) 



Oo 

PC p 



D*°- a ' z+ *4 



Simplify (23) results in, 







Q 0\„ r ^ri^oZ 2 fli S\o +/<T dh.*oZ * z\o 



pc. 



2/7 3/7 



ar 2 /7 s h 



Evaluating this expression gives (24); the buoyancy flux components due to net 
downward temperature flux and bottom heat flux. 



,24) 

In comparing (22) and (24), notice that the depth dependent a causes an 
increased weighting to both the temperature and bottom heat flux components. 
Another difference is the weighting between the two fluxes. The bottom heat 
flux has a larger weight which makes it a larger reduction(addition) in the 
buoyancy flux due to the net downward temperature flux, depending on whether 
AT is positive or negative. This shows that the depth dependent a causes the 
mixed layer to deepen(shallow) at a slightly slower rate for a positive(negative) 
AT than what is found for constant a. 



25 



III. TEST OF PARCEL INSTABILITY THEORY 



A. WATER COLUMN INSTABILITIES 

The response of the mixed layer to the nonlinear effects of a(z) was 
investigated in the previous section, but the possibility of instability below the 
mixed layer needs to be considered. Figure (3) depicts Scott and Killworth's 
(1991) thermistor chain data from a track between Greenland and Iceland. The 
temperature contours show the existence of two chimneys below the surface. 
The chimneys represent events that are clearly below the surface boundary 
layer. Are there features analogous to cumulus cloud formation in the 
atmosphere? Riehl (1979) describes how atmospheric deep convection in the 
tropics requires only a small percentage of the region to be conditionally 
unstable. Is this the same situation seen in deep chimney oceanic convection? 
Scott and Killworth (1991) hypothesized that the chimneys in Figure (3) are 
decoupled from the surface. A chimney could be coupled with the surface if it 
were created by a small instability in a region of conditional instability. A small 
instability from the surface would not be very evident once the chimney event 
started. Thus what looks decoupled after the chimney is developed may not 
show how the surface initiated the event. 

To determine how an instability in a region of conditional instability relates 
to our neutral-parcel instabilities, the model is applied to station data in the 
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Greenland gyre. By displacing a surface parcel downward in the water column, 
the possibility of parcel instability is explored. The surface parcel may be 
cooled and its salinity may also be changed by freezing. The effect on stability 
for such alleged surface parcels is also examined. This procedure 
demonstrates what changes in a surface parcel can cause instability or make 
the parcel neutrally stable/unstable and therefore favorable for a convective 
event. 

B. MODEL DESCRIPTION 

To investigate the possibility of parcel instabilities using in situ data, 
temperature and salinity values from station 44 of MIZEX (Quadfasel et al., 
1988) are applied to the parcel instability model. The temperature and salinity 
data is used in the Millero and Poisson (1981) equation of state to obtain 
density values throughout the water column. The temperature is adjusted for 
the heat of compression with increasing depth by using Bryden’s (1973) 
algorithm to obtain potential temperature (degrees C) as a function of salinity 
(psu), temperature (C) and pressure (bars). 

The in situ density is calculated at each level of the temperature and 
salinity profile. Then moving the surface parcel down the water column profile, 
the temperature is adjusted for compression at each new pressure level using 
the Bryden algorithm. With the compression-adjusted surface temperature and 
surface salinity, the density is calculated for the surface parcel at each pressure 
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level. Each level has an in situ density and a density for the surface parcel 
moved to that level. The change in density is found at each level by 

^ P» = P/iw*v _ Ps**wroor 

When Ap s is negative, the surface parcel would be unstable (more dense) 
compared to its surroundings and would accelerate downward by gravity until it 
reached a depth with equal density or until it encountered the bottom. 

Station 44 will be examined to determine if the situation described above 
ever occurs. The results are seen in Figure (7), a plot of depth versus Ap. 

The solid line shows the non-compression adjusted temperature values for the 
change in density at each level as depth increases. The dashed line shows the 
adiabatic temperature adjusted surface parcel changes in density as depth 
increases. The profiles do not exhibit any parcel instabilities for the displaced 
surface parcel since the values of Ap never become negative. Figure (7) also 
shows that the temperature adjustment is a stabilizing effect for the profile; the 
adiabatic temperature corrected Ap curve is more positive through the water 
column. 

1. SENSITIVITY TO DEPTH 

The procedure described is then modified to begin with a parcel at 
the bottom of the mixed layer vice a surface parcel. No parcel instabilities were 
noted. Other levels are also examined, and no parcel instabilities were 
discovered. Figure (8) is an example of this procedure. The dot-dash line 
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Figure 7. Change in density curves for Station 44: The dotted curve 
represents adiabatic temperature adjustment. The solid curve is not adiabatic 
temperature adjusted. 
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CHANGE IN DENSITY (KILOGRAMS PER CUBIC METER) 




Figure 8. Change in density curve: Testing for depth sensitivity, solid being 
non-adiabatic, dotted being adiabaticly adjusted and dot-dash is a 500 m parcel 
moved in the water column. 
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shows a parcel from 500 meters depth moved to the surface and then moved 
down the water column. Below 500 meters the 500 meter parcel never became 
unstable. 

It is worth noting that Gascard (1991) stated that chimney formation 
only occurs were the overall stratification can be destabilized and destroyed 
from the interior. He further stated that vertical mixing by turbulent entrainment 
triggered from the surface is ineffective. Figure (8) shows an example of a test 
for destabilization of the stratification from the interior, i.e. 500 meters in this 
example. However, looking at the above parcel instability model for the data in 
the Greenland Sea given by Gascard (1991) shows no parcel instabilities using 
water from intermediate levels. Figure (9) demonstrates how a 1000 meter 
parcel moved in the water column did not cause instability either. 

2. SENSITIVITY TO SURFACE TEMPERATURE 

The next test was to assume that the surface was cooled in steps 
down to the freezing temperature by either heat loss to the atmosphere or by 
advection. By using Millero's (1978) equation for the freezing point of 
seawater, the maximum amount of cooling that could be imposed on the 
surface is calculated. At station 44, the surface temperature was already - 
1.797 C and could only be cooled by 0.1 C to reach the freezing point. 

The first step is to apply the parcel instability model with station 44's 
original surface temperature. Each consecutive test is applied to the same 
station 44 water column profile but with the surface parcel temperature 
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CHANGE IN DENSITY (IN SITU SOLID, PRESS. ADJ. DOT, 1000M DOT/DASH) 
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Figure 9. Change in density: Test of Gascard’s Greenland Sea data for 
internally caused instabilities. 
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CHANGE IN DENSITY (KILOGRAMS PER CUBIC METER) 




Figure 10. Contour of Ap(0,Z): Sensitivity to reducing surface temperature to 
freezing. 
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temperature in C 



reduced by 0.02 C at the start of each iteration until the freezing temperature 
was reached. Figure (10) shows the contours of Ap(0,z). Viewing Figure (10) 
from right to left shows the result of decreasing station 44‘s surface parcel 
temperature to the freezing point. 

Figure (1 1) shows an expanded view of the upper 750 m for the 
same contour intervals. The only instabilities demonstrated are in the upper 50 
m of the water column, the mixed layer. One interesting feature to note is a 
zone of maximum stability between 500 m and 1500 m that is strongest at the 
original surface parcel temperature and becomes less dominant as surface 
parcel temperature approaches freezing. This can is evident from the 0.08 Ap 
contour in Figure (11). 

3. SENSITIVITY TO SURFACE SALINITY 

With no parcel instabilities apparent in the water column when the 
temperature was reduced to the freezing point, the effects of surface salinity 
changes are examined. An increase in salinity is associated with freezing after 
the surface water is cooled to the freezing point. Then the change in salinity 
(6S) times the thickness of the ice (h) is proportional to the change in the ice 
thickness (6h) times the difference in surface salinity (S) and ice 
salinity (S,), 

bS*h*bh*(S-S) 

The salinity of the ice is approximated with Maykut's (1985) ice salinity data. 
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Figure 11. Contour Plot of Ap(0,Z): Expanded view of the upper 750 meters of 
sensitivity to reducing surface temperature to freezing. 
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temperature in C 
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Figure 12. Contours of Ap(S,z): Sensitivity to increasing salinity by .02 psu 
with the surface temperature at freezing. 
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salinity in PSU 



An ice salinity value of 12 psu for first year ice of initial thickness at 5 
cm or less will give approximately 0.1 psu increase in salinity for 1 cm of ice 
growth at the surface of the water column. This does not account for the 
salinity increases due to mixed layer turbulent mixing or advection of colder 
more saline water into the surface water. Anderson's (1961) data for ice growth 
rate gives a 10 cm/day ice growth rate for a storm with air temperatures at -40 
C. Thus an increase of 0.1 psu for the surface water is reasonable to use in 
the parcel instability model. 

Solving the parcel instability model for the surface parcel at freezing 
temperature and increasing the salinity by 0.1 psu in 0.02 psu increments gives 
the contour plot of figure (12). The graph is salinity vs. depth, plotting the 
change in density between the surface parcel moved to depth as salinity 
increases by 0.02 psu increments, Ap(salinity, depth). As the salinity is 
increased from its original value at the surface of 34.74 to 34.80 psu and 
greater, the parcel instability (negative values of Ap) are predicted. Also, values 
greater than 34.84 psu show parcel instability to the bottom of the water 
column. The second to the last iteration (second to last contour on the right) 
shows the neutral stability curve for a parcel moved from the surface to depth 
without any change in the density between the surface parcel at the freezing 
point and the in situ water density. 

Figure (13) shows an expanded view of the upper 750 m for the 
same contour intervals. This shows how the mixed layer instability grows very 
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Figure 13. Contour Plot of Ap(S,z): Expanded view of upper 750 meters, 
sensitivity to increasing salinity with surface temperature at freezing. 
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salinity in PSU 



quickly from the onset of the salinity increase. Examining the values of 
instability near the salinity value of 34.76 psu shows at what value parcels 
should "break out" of the mixed layer in deep convection events. Thus the 
model demonstrates that a parcel instability can occur for the station 44 profile 
provided that the surface is cooled to freezing and as little as 1 cm of ice forms. 
The salinity injection and surface temperature changes at the surface can 
create a parcel instability that can cause convection to the bottom of the water 
column. 
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IV. CLIMATOLOGY 



A. CLIMATOLOGY AND PARCEL INSTABILITY 

Is there a region in the Greenland sea that can be considered preferential 
for parcel instability events leading to chimneys? With the assumption that 
parcel instabilities depend on a depth dependent thermal expansion coefficient 
the Levitus(1982) Climatological Atlas of the World Oceans is used to 
investigate the Greenland Sea. The Levitus (1982) data set contains 
temperature and salinity at standard oceanographic observation levels on a one 
degree latitude and longitude grid for the world oceans. The data are averages 
of objectively analyzed station data, mechanical bathythermograph and 
expandable bathythermograph between the surface and a maximum depth of 
5500 meters. Annual and seasonal data were provide, but only the winter 
season data is used for the parcel instability analysis here. 

B. CLIMATOLOGICAL MODEL TEST 

Again using the depth dependent thermal expansion coefficient and the 
principles of the parcel instability model, the temperature and salinity values 
from 55 N to 85 N in latitude and 5 E to 20 W in longitude will be used to 
determine the stability of the region. The surface values of temperature and 
salinity from the data set do not represent a very realistic mixed layer due to 
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the averaging process, but it is assumed they provide an adequate first guess 
field. A depth of 500 m will be the comparison depth for the instability test. 

This depth is used to help cover the largest area of the data set without being 
bottom limited. The 500 m depth is also a reasonable chimney depth evident in 
various data sets. 

The model will use (25) to determine the buoyancy jump (Ab) as a 
function of the 500 m depth to create a grid of the change in buoyancy ( Ab) 

< 25 > 

values to contour. The values subscripted with 500 indicate the value at 500 m 
and the subscript surf indicates the surface value. Again, in this test (3 will be 
kept constant. 

Figure (14) represents the Ab values plotted on the latitude and longitude 
grid. The area encircled by the 0.1 cm/s 2 contour represents an area close to 
neutral stability. Comparison of this area to the charts in Figures (1) and (2) 
show that this is highly correlated to the central region of the Greenland Gyre. 
The areas on Figure (14) without a contour value represent land or the water 
depth was too shallow for comparison with the 500 m depth. 
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Figure 14. Contoured Plot of Change in Buoyancy: Contours of the Ab 
between the surface values and a depth of 500 meters. 
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LONGITUDE (WEST IS NEG.) 



This instability test supports the theory that the Greenland Gyre is the 
most likely candidate for chimney formation in the Greenland Sea region. It 
helps provide a bounds for further data collection and observations of chimney 
events. 
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V. SUMMARY 



This study has (1) investigated the dynamics of a neutral-parcel model 
with a depth varying thermal expansion coefficient (a), and (2) tested for parcel 
instability in the Greenland Sea, using both actual profiles and climatological 
data. 

In the neutral-parcel model, the effects on mixed layer dynamics for forced 
and free convection and the combination of both were analyzed for what added 
dynamical effect the depth dependent a(z) had on the mixed layer in each test. 
The depth dependent a(z) changes significantly mixed layer deepening. For all 
three cases of convection, depth dependent a(z) caused the mixed layer to 
deepen(shallow) at a slower rate given a positive(negative) AT than what was 
found for a constant a. 

The parcel instability model using a(z) and adiabatic temperature 
increases due to parcel movement down the water column was tested for 
sensitivity to parcels displaced from interior depths, sensitivity to reducing the 
surface temperature to freezing and sensitivity to increases in the surface 
salinity. Gascard's (1991) destabilization from the interior could not be 
supported using the parcel instability model, demonstrated by the test for depth 
sensitivity using data from Gascard (1991). Reducing the surface temperature 
to freezing showed the onset of freezing reduced the stability and brought it 
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closer to neutral stability. Cooling surface temperature to freezing did not by 
itself lead to a state of instability for most of the profiles tested. When the 
surface was at freezing and enough ice formation occurred, then the increase 
in salinity resulted in the possible onset of parcel instability. 

Using the depth dependent thermal expansion coefficient and adiabatic 
temperature increases, an assessment of the possible instability in the 
climatological data of the Greenland Sea was made. The area of least stability 
for this test was the center of the Greenland Gyre. 
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